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The main goal of this paper is to extend and apply the principle of maximum entropy (Max- 
Ent) to incomplete quantum process estimation tasks. We will define a so-called process entropy 
function being the von Neumann entropy of the state associated with the quantum process via 
Choi-Jamiolkowski isomorphism. It will be shown that an arbitrary process estimation experiment 
can be reformulated in a unified framework and MaxEnt principle can be consistently exploited. We 
will argue that the suggested choice for the process entropy satisfies natural list of properties and it 
reduces to the state MaxEnt principle, if applied to preparator devices. 
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I. INTRODUCTION 

Physical objects and processes are described by pa- 
rameters that are directly, or indirectly, accessible exper- 
imentally and represent the maximal knowledge about 
physical systems (according to physical theory used). In 
quantum theory (see for instance [1-3]) the complete in- 
formation (knowledge) is represented by the concepts of 
quantum state (normalized positive operator), quantum 
observable (normalized positive operator valued mea- 
sure) and quantum channel (completely positive linear 
trace-preserving map). One of the main characteristics 
of quantum system is its dimension d, i.e. the maximal 
number of mutually perfectly distinguishable states (in 
a single run of the experiment). These states form an 
orthogonal basis of the associated complex Hilbert space 

n. 

An arbitrary quantum state is described as a posi- 
tive trace-class linear operator with unit trace acting on 
the Hilbert space, g : H — » H : g > 0,Trg = 1, i.e. 
a density operator, or a density matrix. The number 
of independent real parameters determining the quan- 
tum states scales as A4tate = d? — 1. The quantum 
processes/operations correspond to completely positive 
trace-preserving linear maps defined on the set of all 
linear operators including the set of all states S{TL). 
The number of independent real parameters determin- 
ing the particular quantum operation equals A/process = 
d 2 (d 2 — 1). Quantum measurements give us probabil- 
ity distributions over the set of all possible outcomes 
{xi, . . . ,xl}, where L is some positive integer. In the- 
ory, the measured probabilities pj are determined by the 
Born's rule pj = TrgFj, where Fj is a positive operator 
(quantum effect) corresponding to outcome Xj . These op- 
erators form the so-called positive operator valued mea- 
sure (POVM), i.e. Fj are positive (Fj > 0) and they sum 
up to identity operator (Ylj Fj = I). The number of pa- 
rameters specifying POVM depends on the total number 

Of Outcomes L and equals Admeasurement — {L — l)d 2 . 

The goal of quantum tomography is to estimate and 
fix all these parameters [4-6]. However, already for small 
systems (in dimension) the number of parameters is in- 
creasing rapidly, especially for quantum channels [7—11]. 



It seems that the complete knowledge about quantum 
objects is not a very realistic dream and experimentally 
we will not be able to perform all the desired tests [12- 
14]. Fortunately, there are situations in which even the 
knowledge of only few parameters enables us to make 
reasonable and nontrivial predictions about the behavior 
and properties of the system. A typical (classical) exam- 
ple is the equilibrium thermodynamics in which only few 
parameters are used to describe the complex behavior of 
a system of approximately 10 23 degrees of freedom. Our 
aim is to describe the properties of quantum objects as 
honestly as possible even in cases when the complete in- 
formation is not available. In particular, in this paper we 
will focus on incomplete quantum process tomography. 

In Section II we will define the concept of process mea- 
surement and shortly describe the idea of quantum pro- 
cess tomography. The maximum entropy principle is de- 
scribed in Section III and also the idea is extended to 
process estimation problems by introducing the concept 
of process entropy. Finally, in Section IV the MaxEnt 
procedure is applied to particular examples of incomplete 
ancilla-free estimation of qubit channels. 



II. QUANTUM PROCESS MEASUREMENT 

A general quantum process tomography experiment 
consists of a test state g that is transformed in some 
specific procedure V involving the unknown channel £ 
into a state g' , and a measurement (POVM) M. per- 
formed on the state g' . This framework includes all the 
possible strategies [5, 6, 14] via which the parameters 
of quantum channels £ are accessible. We will define a 
process measurement as a particular choice of the test 
state g, of the procedure V and of the measurement M.. 
Generally, the procedure V is composed of an applica- 
tion of some known quantum channels on the test state 
and one usage of the unknown channel £ acting on d- 
dimensional quantum system (qudit). That is V itself 
is a quantum channel that can be written as a product 
V = V in o (F &ac ®£ )o r P out , where V m ,V OVl x can be under- 
stood as being parts of the preparation of the initial state 
g, and of the final measurement A4, respectively. Con- 
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sequently, without loss of generality we can assume that 
V = T anc ® £, where is a known quantum channel 
acting on some ancillary system and also can be included 
as being a part of either preparation of g, or measurement 
performed. A process measurement is called ancilla-frce 
if either the initial state g is factorized, or the ancillary 
system is trivial. Otherwise the process measurement is 
ancilla-assisted and g' = l anc <g> £[g}. 

For example, consider £ is a qudit quantum chan- 
nel and the test state g = is a maximally entan- 
gled state of two qudits (*+ = g^liX^I ® 
The unknown channel is applied only on second of the 
qudits while the first one is transformed trivially, i.e. 
g' = u)£ = 1 (g> £ [6]. Performing the information- 
ally complete POVM (resulting in complete specification 
of cog) the channel £ can be uniquely identified, because 
the mapping £ i— > uj£ — J\£\ is the well known Choi- 
Jamiolkowski isomorphism [15, 16] between the set of 
quantum qudit channels and set of quantum states of two 
qudits. Hence, via general POVM measurements of the 
output state we can acquire either complete, or partial, 
knowledge on the channel. Let us note that individual 
ancilla-free process measurements cannot be information- 
ally complete, but they can be combined together to gain 
the complete information. In the following sections we 
will concentrate onto situations in which the collection 
of process measurements provides us with partial infor- 
mation, only. 



III. PRINCIPLE OF MAXIMUM ENTROPY 

Maximum entropy (MaxEnt) principle was originally 
introduced in statistics in order to estimate a probability 
distribution providing that only partial information on 
that probability is available [17]. There are many proba- 
bility distributions compatible with the given constraints 
and our aim is to choose one of them that in some sense 
represents our knowledge the most honestly. This choice 
cannot be logically derived and some additional principle 
must be introduced. Using the results of information the- 
ory [18] on the uniqueness of Shannon entropy, one can 
argue that [17, 19] the probability distribution maximiz- 
ing the Shannon entropy is the best choice we can make. 
Such probability maximizes the uncertainty (measured 
by entropy) and, intuitively, also our predictions about 
the unspecified parameters are as uncertain as possible. 
That is, a conclusion based on MaxEnt principle is in- 
troducing as little additional information as it is possible 
[19]- ^ 

This idea was generalized to the domain of quantum 
state tomography [20, 21] by using the concept of von 
Neumann entropy [1] S(g) = — Trglogg, which is consid- 
ered to be the quantum extension of Shannon entropy. A 
state observation level 0„ is defined as a set of n (n < d 2 ) 
mean values {/i, . . . , /„} of linearly independent opera- 
tors {F\ , . . . , F n } related to unknown state g via the trace 
rule fj — TrgFj = (F} g . If the operators {iq, . . . , F n } 



form a POVM the numbers fj represent the measured 
probabilities. The answer to incomplete state tomog- 
raphy problem based on maximum entropy principle is 
given by the following equation 

g = argmax{S(g)|£ e S(H), fj = TrgFj, j = 1, . . . ,n} . 

The following state is the formal solution of the MaxEnt 
problem [17, 21] 

e=|exp(-E A ^)' ^ 

j 

where Z = Tr[exp(— ^\ XjFj)] and Aj are Lagrange mul- 
tipliers fixed by the system of equations 

fj = TrgFj = -— InZ(Ai, . . . , A„) . (3.2) 

For example, consider a two-level quantum system 
(qubit) and observation levels 

Oi ={<**>„} 

C 2 = {{Oy) e ,{(Tz) e } 

C 03 = {(Crx) e , {<Ty) ei {(T z ) e } . 

A qubit state can be expressed in a so-called Bloch sphere 
picture as g = ^(1 + r- a) with Tj = Tr go~j. The MaxEnt 
principle applied for 01,02 sets mean values Tj of all 
the unobserved operators to zero. That is, for 0\ = {z} 
we get g = \{I + z<r z ) and for Oi — {y, z} the MaxEnt 
estimation gives g = \{I + yoy + za z ) . Observation level 
03 provides complete information about the quantum 
state, hence the principle of maximum entropy is not 
needed in this case. 

Clearly, there is a problem if we consider similar in- 
complete estimation task for processes, namely, which 
entropy should be maximized? Unlike quantum states 
the quantum channels are lacking some concept of en- 
tropy, or uncertainty. In fact, what does it mean that a 
quantum process is uncertain? Our goal is to introduce 
a suitable concept of a channel/process entropy S proc (£) 
and investigate its properties. Before analyzing different 
choices let us discuss some (intuitive) properties of the 
process entropy. 

1. Uncertainty of unitary channels. Without any 
doubts the unitary channels play a very specific 
role among all quantum processes. For unitary pro- 
cesses the interaction of the system with its envi- 
ronment is trivial. The physical invertibility is the 
unique and characteristic property of the unitary 
channels. In some sense the channel entropy should 
reflect how much noise the channel introduces. Uni- 
tary channels are noiseless and in what follows we 
will assume that the channel entropy is invariant 
under unitary preprocessing (V) and unitary post- 
processing (14), i.e., Sp r oc(£) = S pioc (U o £ o V). 
It follows that all unitary channels have the same 
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value of uncertainty that can be set to zero. More- 
over, we do require that S proc (£) = implies that 
£ is unitary. 

2. Uniqueness of maximum. For the purposes of in- 
complete process estimation exploiting the Max- 
Ent principle it is necessary that the maximum is 
unique. Hence there must be a unique channel, for 
which the uncertainty is maximal. This channel 
should be the result of the incomplete estimation 
if no data are available, i.e., when the observation 
level is trivial, O = 0- Because of the unitary in- 
variance the channel must be invariant under uni- 
tary preprocessing and postprocesing, i.e., £ max = 
U o £ max o V. Only the channel mapping the whole 
state space into a total mixture (g i/) is invari- 
ant in this sense. It is argued in [22] that this chan- 
nel is indeed the average channel over all possible 
qubit channels. To guarantee that in any process 
measurement the maximum is unique it is sufficient 
that the process entropy is a concave function, i.e., 
5proc(A£i + [l— A]£ 2 ) > XS pTOC (£ i)+[l— \]S pmc (£2)- 

3. Universality. This is not a condition on the con- 
cept of process entropy itself, but rather on the 
general possibility to employ such principle once 
we agree on a suitable measure of channel entropy. 
It is important that the maximum entropy princi- 
ple is applicable for all process measurements. We 
shall discuss this issue later in more details. 

In summary, a channel entropy is some concave func- 
tion (defined on the set of quantum channels) achieving 
its maximum for the complete contraction to the total 
mixture and vanishing only for unitary channels. A nat- 
ural choice of process entropy seems to be related to the 
concepts of quantum channel capacity [23-26] . Quantum 
capacity quantifies the degree of preservation of quantum 
states during the transmission and this value is different 
for different unitary transformations. On the other hand, 
the classical capacity is maximal also for noisy channels. 
For example, phase-damping channels g — > diag[g>] max- 
imize the transmission of classical information over the 
quantum channels. Because of these properties, the ca- 
pacities are not appropriate candidates for the definition 
of process entropy usable in incomplete process tomog- 
raphy tasks. 

A. Choi-Jamiolkowski process entropy 

The Choi-Jamiolkowski isomorphism provides us nat- 
urally with a notion of channel entropy. It uniquely asso- 
ciates a quantum state cog = (Z"(g> £)[*+] with a quantum 
channel £, hence we can adopt the von Neumann entropy 
of ui£ as being the channel entropy of £ [27, 28]. Con- 
sider a quantum channel on rf-dimcnsional system (qu- 
dit). Providing that for each process measurement we 
are able to define uniquely a state observation level O n 



given by mean values Xj = TrXjUg of n linearly inde- 
pendent Hcrmitian operators Xj, the MaxEnt problem 
for processes can be formalized as follows 

£ — argmaxS^wg) 

where the maximum of von Neumann entropy S{u>s) 
is taken over all states uj£ E S(H <g> Ti) satisfying the 
constraints Tr2a>£ = 4/ and Xj = TruisXj for all 
Xi, . . . , X n e O n . The resulting state uj£ determines the 
quantum operation £ uniquely via the inverse relation 

£[g] = dTr anc [{ g T ® I)u> £ ] . (3.3) 

In what follows we will investigate the process entropy 
given as the von Neumann entropy of the state tog asso- 
ciated with the channel £ via Choi-Jamiolkowski formal- 
ism. 

It is straightforward to see that only for unitary chan- 
nels the states tog are pure and hence S proc (£ ) = S(u>£) = 
only for unitary channels, £ — U. Moreover, uni- 
tary channels do not change the entropy of u)£, i.e., 
Situs) = S(ujuoSov)- The concavity of S ploc follows from 
the concavity of von Neumann entropy and the maxi- 
mum is achieved for uj£ = -^1 that is associated with 
the channel mapping the whole state space into the max- 
imally mixed state, £ : g 1— > ij. In summary, the process 
entropy 

S woc {£) = -Trw f logw £ (3.4) 

satisfies all the desired properties we have discussed pre- 
viously. The only open issue is its applicability in general 
process measurement. 

Choi-Jamiolkowski isomorphism is associated with a 
specific process measurement using as the test state a 
maximally entangled state of two qudits. Second qudit 
is sent through the unknown channel while the first one 
is evolving trivially to obtain the state X ® £ = , 
that is estimated in some state measurement described 
by POVM. In this case the process observation level can 
be defined as the following set of mean values 

Or c = {xi,...,x n }, (3.5) 

where 

x i = Tr -F>£ = (Fj)T®£[n, +] = (Fj) £ . (3.6) 

Because of the identity Tr2Wg = 2^ ^ ne process observa- 
tion level is equivalent to a state observation level 

Cn+d 2 -l = {xi, ■ ■ ■ ,x n ,0, . . . ,0} . 

The added zeros represent the mean values of d 2 — I op- 
erators / <B> Aj , where Aj are tracelcss Hermitian qudit 
operators forming a basis of the set of traceless Hermi- 
tian qudit operators, i.e., the general qudit state can be 
written as g = g (/ + r ■ A) . To be more precise we assume 
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that the operators F\ , . . . , F n are linearly independent of 
operators I ® A l7 . . . , I ® A d 2_ 1 . 

What if the test state is not the maximally entan- 
gled one? Is it possible to interpret the measured values 
as linear constraints on the state Sl £ defined by Choi- 
Jamiolkowski isomorphism? Let us note that the lin- 
earity is crucial, because we implicitly assume that the 
constraints representing the incomplete information are 
linear, which guarantees that the set of possible solutions 
is convex and, hence, the entropy has a unique maximum. 



B. General quantum process experiment vs 
Choi-Jamiolkowski isomorphism 

Consider a general test state SI of the qudit and an 
arbitrary ancilla system. We will show that there exist 
a completely positive linear map An : B(Hd) — * B(Ti anc ) 
such that An ® ^[*+] = fi. A general pure state 

I*) = Eaj $ a>)anc ® U) {a = 1, ■ ■ ■ , D; j = 1, . . . d) 

can be written as |$) = A$ ® where the operator 

A$ : — > 7Y anc is defined as A$ = VdJ2 a .j ;®aj\a)(j\- 
A general mixed state SI can be written as convex com- 
bination of pure states SI = ^ fc \k\&k)(®k\, hence SI — 
J2 k A fc ( <gi /) * + (At (gi J) = (An® !)[*+]. Since the val- 
ues Afc are positive the transformation An is completely 
positive. Moreover, for each state S} the linear map An is 
unique. Hence, for a general test state SI the mean value 
of an Hcrmitian operator F can be expressed as follows 

(-4s2®£)[ > I'+] 

= ((Ah®l)[F]) {m£)[ * +] , (3.7) 

where Aq is a dual mapping to An (Hciscnberg picture). 
As a result we get that an arbitrary ancilla-assisted pro- 
cess measurement can be rewritten within the framework 
of process measurement using the maximally entangled 
test state and measuring a suitable Hermitian opera- 
tor Aq®T[F], hence, the maximum entropy principle de- 
fined via Choi-Jamiolkowski entropy can be consistently 
employed in all incomplete process measurements. 

In what follows we shall analyze the ancilla-free pro- 
cess measurement, hence only the qudit itself is used to 
probe the action of the quantum channel. In fact, this 
can be considered as an ancilla-assisted problem with 
a factorized test state SI = £ anc ® Q, and a factorized 
measurement resulting in the mean value of the operator 
of the form J anc ® F. Consider is a pure factorized 
state |$) = |(ys) a nc ® | -0) - Then the operator A$ takes 
the following form A$ = \(p)(ip*\, where \ip*} is a com- 
plex conjugated state, i.e., (k\ip*) = (k\ip) for all basis 
vectors in which the maximally entangled state \& + 
is defined. It follows that for general factorized state 
^ = E a ,i ^aWlv>a)v>a| ® \^l){i>l\ tnc transformation An 
is expressed via Kraus operators A al — y/dX a ni\ifi a )(ip*\- 
Therefore, according to Eq. (3.7) the ancilla-free pro- 
cess measurement of 7 anc ® F can be considered as an 



ancilla process measurement with the maximally entan- 
gled state * + being the test state and a measurement of 
X = (A n ® 1) [/anc ® F] , i.e., 

x = ^(A^®/)(/ anc ®F)(A aji ®/) 

a,l 

= dJ2 X a»l\4>t)(<P a \<P*)(llj!\®F 
a.l 

= d(£»im(m®F 
i 

= dg T ®F, 

where we used that g = X^mI^X^Ii Cane = 
J2 a ^a\fa)(fa\ an d Q T is the transposed matrix g with 
respect to basis {\k)}. 

We have shown that measuring the outcome associ- 
ated with F in the ancilla-free process measurement is 
equivalent to measuring dg T ® F in the process mea- 
surement with maximally entangled state \&+, where g 
is the ancilla-free test state. It means that the ancilla- 
free process observation level consisting of mean val- 
ues (F 1 ) ei . . . , (F n ) Bn is equivalent to 0P roc = {(dgj ® 
Fi)s, . . . , (dgT ® F n )s}. 

C. States as preparation channels 

Preparation devices play a completely different role 
than quantum channels. However, formally, they can 
be understood as mappings that transform an arbitrary 
input state into a fixed output state £. In this sense 
preparation channels £^ form a very specific convex sub- 
set of quantum channels. Let us apply the proposed max- 
imum entropy based process tomography to preparation 
channels, i.e., to preparators. The process measurement 
is ancilla-free consisting of all linearly independent test 
states Qj (j = 1, . . . , d ). and measurement of the mean 
value of Hcrmitian operator F. According to previous 
paragraph the process observation level is described as 
O p d T c = {(dgj ® F) £( , . . . , (dg T d2 ® F) £ J. The Choi- 
Jamiolkowski entropy of the channel £^ equals (up to a 
constant) to the von Neumann entropy of the state £, be- 
cause Z®£ e [4> + ] = 3-T®£ implies S(w £( ) = log 2 d+S(£). 
Moreover, because of the identity 

(dgj ® *W 4[ * +] - (dgj ® F)i m = (F)t , (3.8) 

it follows that finding a channel with the maximal Choi- 
Jamiolkowski entropy is equivalent to finding a state 
maximizing the von Neumann entropy. As a result we 
get that the process MaxEnt procedure, if applied to 
preparators, reduces to the state MaxEnt procedure. 
That is, the proposed Choi-Jamiolkowski process entropy 
is a consistent extension of the von Neumann entropy. In 
particular, the MaxEnt principle for states can be con- 
sidered as being a special case of the MaxEnt principle 
for channels. 
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IV. EXAMPLES 

In this section we shall present few examples of incom- 
plete quantum process estimation for ancilla-free process 
measurements of a qubit channel. 

A. Ol IOC = {{2q t ®<j z ) £ } 

In this case the collected data provides us about in- 
formation on the mean value of an observable a Zl hence 
the experiment gives us single value to = (o , z }s[g\- Unfor- 
tunately, even in this simplest case we cannot give (see 
Appendix A) an analytic solution in its whole generality. 
In particular, we found the solutions in following cases 

g=\l : £[Z\ = \(I + m<T x ), ( 

where g = \ (I + f ■ a) and £ = 5 (7 + t ■ a) . It is inter- 
esting that for pure test state the estimated channel is 
not unital even if there are unital channels satisfying the 
constraints. 

An alternative method for incomplete process estima- 
tion was described in [22]. It is based on a different ad 
hoc rule demanding that no additional information about 
unobserved measurements (those completing the incom- 
plete process observational level) is introduced. In par- 
ticular, for states 77 orthogonal (in Hilbcrt-Schmidt sense, 
i.e., Tvqg = 0) to given test states, the mean values are 
completely random, i.e., they are transformed into the 
total mixture (77 — > ij). Hence the entropy of output 
states for unmeasured inputs is maximal. It means that 
if possible (meaning there is no contradiction with the 
data, or theory) the total mixture is preserved. Other- 
wise an optimization procedure minimizing the average 
distance from the total mixture is needed. This method 
was analyzed only for qubit channels and for ancilla-free 
process measurements. The extension of the method to 
all process measurements will require introduction of ad- 
ditional rules. Let us compare the method proposed in 
[22] and the one proposed in this paper. 

If measuring o~ z and finding m = ±1 the output state 
must be pure and it corresponds to an eigenvalue of a z . 
In both mentioned scenarios we know the solution pro- 
viding our knowledge consists of complete information of 
the action of the channel on the pure test state, thus, we 
know that £ : \ip) 1— > \±z), respectively. As it was argued 
in the work [22] the estimated transformation should map 
the whole Bloch sphere into the line connecting north and 
south pole, i.e., t — > r = (0,0, ±t 2 ). However, the pro- 
posed MaxEnt estimation procedure gives different re- 
sult. In particular, t -> F = (0, 0, ±(1 + t z )/2). This 
transformation is not unital, but the total mixture is 
mapped to the state 7? = (0,0,1/2). A state t = —r 
orthogonal to the test state r is transformed as follows 

Scheme in [22] : to = 1 £ cst : —r^> —r 
MaxEnt: m=l L t : -?n ^ ' ' 



As we see in this case both methods transform orthogonal 
(in Hilbcrt-Schmidt sense) states to \ip) into the total 
mixture, but for MaxEnt procedure also the orthogonal 
(in Hilbert space sense) state \tp±) is mapped into the 
total mixture. In our opinion this feature (except the 
universality) justifies the usage of MaxEnt procedure in 
comparison with the scheme described in [22]. In fact, 
the uncertainty introduced by the estimation procedure 
on perfectly distinguishable (orthogonal) states from the 
test states should be as maximal as possible. And this is 
not the case for the method used in [22], for which the 
estimated channel preserves the orthogonal state. 

B. Ol TOC = {</ ® a x } £ , (I ® o v )s, (I ® a z ) £ } 

Consider a situation that the unknown qubit channel 
is tested by the total mixture and the complete tomogra- 
phy of the output state is performed, i.e., mean values of 
a Xl (j yi (j z are known. The corresponding state observa- 
tion level is Oq = {{<j®I) e , {I®a) e } = {0, to}, for which 
the solution is presented in Appendix B. In such case the 
proposed MaxEnt process tomography procedure leads 
us to the channel 

£ cs t ■ £"-► Qa = ^(I + m-a), (4.3) 

hence, the whole Bloch sphere is contracted into a single 
point to. As a result we get that if the total mixture 
is used to probe the channel action then according to 
MaxEnt procedure all the states are mapped into the 
output state go — £[^I]- In this case both the discussed 
procedures are giving the same estimation. 

C. Of oc = {(/ ® a z ) £ , (2(\x)(x\) T ® a z ) £ , (2{\y)(y\) T ® 
a z ) £ ,(2{\z)(z\) T ®a z ) £ } 

In this case the process is probed with four test states 
(total mixture and positive eigenvectors of <j x ,a y , o~ z 
forming a vector of pure states if), but only zth compo- 
nent of the Bloch vector of the output state is measured. 
07 = {</ ® <Tz) e , (2rf ® <r z ) e , (a ® 7)J = {z, (, 0} is the 
corresponding state estimation problem and z, £ are the 
experimentally identified mean values. Information en- 
coded in these parameters can be equivalcntly rewritten 
into the form O7 = {{I ® & z ) e , (<? ® cr z)g, (<? ® I) e ] = 
{z,C',0}, where £j = Q - z. 

It is shown in [21] that for such state observation level 
the estimated density matrix reads 

uj = j (i ® I + zl O a z + (f • a) <g> a^j . (4.4) 

Hence, the process is described by the following state 
transformation (g — > £ cst [g\) 

f-» ? = (0, 0, z + C' • t) , (4.5) 
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i.e., g' = j[I + {z + Cf ■tycjz]. As in all previous cases, also 
in this case the whole state space is mapped onto a subset 
of the line connecting states | + z) and | — z). However, 
in this case the final state depends also on parameters 

V. CONCLUSION AND DISCUSSION 

We have addressed the problem of incomplete process 
estimation based on maximum entropy principle [19]. In 
general the maximum entropy principle is an intuition- 
based ad hoc principle related to quantification of igno- 
rance contained in probability distributions that seems to 
agree with our experience. This ignorance measured in 
entropy can be extended to domain of quantum states by 
introducing the von Neumann entropy. Our attempt here 
was to develop similar approach for processes. We argued 
that capacities are not good candidates for quantifying 
the uncertainty of quantum channels and we exploited 
the Choi-Jamiolkowski process entropy defined as 

S P roc(£) = -Tr[u £ logw £ ] , lo £ = (J® £)[*+] , (5.1) 

where ^+ is the maximally entangled state. In particu- 
lar, we showed that the suggested concept can be univer- 
sally applied in all possible process measurements. The 
procedure is demonstrated on three incomplete ancilla- 
free estimation problems of a qubit channel: i) pure test 
state and projective measurement, ii) the total mixture 
as the test state and complete tomography of the output 
state, and iii) four test states and the same projective 
measurement. 

We have shown that unlike the concepts of capacity of 
quantum channels the process entropy defined above is 
compatible with the following properties: 

1. Uniqueness of maximum: unique maximum for the 
channel contracting whole state space into the total 
mixture. 

2. Minimum: minimum is achieved only for unitary 
processes. 

3. Unitary invariance: invariant under unitary trans- 
formations, i.e., S pTOC (U o £ o V) = S proc (£) for all 
unitary transformations U, V. 

4. Concavity: function S proc (£) is concave, i.e., 
S proc (p£ + qT) > qS pioc (£) + pS pi0 c(^)- 

The proposed Choi-Jamiolkowski process entropy serves 
as a very valuable tool in incomplete process tomography 
deserving future testing and investigation. Moreover, as 
it is shown in Section III.B, the proposed process entropy 
principle, if applied to state preparator devices, is equiva- 
lent to the state entropy principle based on von Neumann 
entropy. The key feature discussed in this manuscript is 
the universality of the proposed procedure following from 
the unification of all process measurements described in 



Section III.B. This idea goes beyond the applications in 
incomplete process estimation and is further developed 
in [30]. 

Recently, Olivares et al. [29] proposed and analyzed 
a state estimation problem combining incomplete in- 
formation with some nontrivial apriori knowledge. In 
their approach the maximization of entropy is replaced 
by minimization of Kullback relative entropy S(q\qo) = 
Tr[g(log g — log go)] with a bias go representing the prior 
knowledge. This approach can be directly extended 
to the case of channels by introducing the quantity 
S , (a>£||cjo) = — Tr[aj,s(log u>£ — log luq)} with £ o playing the 
role of prior information. If we set £q to be the state space 
contraction into the total mixture (i.e. ujo = ^r-Oi th cn 
Tr[u>£ logcjg] = —logd 2 Tru)£ = — 21ogd. Consequently, 
S^c^H^o) = 21ogd — S(uj£) and the biased estimation 
problem reduces to the unbiased maximum process en- 
tropy estimation. 

Let us give a simple example based on the observa- 
tion level discussed in Section IV. A. Suppose that out 
of the performed measurement we acquire the informa- 
tion |0) i ^ |0). We shall consider three different pri- 
ors: i) identity channel £ = X; ii) diagonalisation chan- 
nel £o = diag transforming each state into its diagonal 
form in the basis |0),|1); iii) £o = \{T +U X ), where 
M X [Q = o~x£,o~ x . Let us note that 5 , (g||go) is finite only 
if the support of q is included in the support of qq. For 
the case of identity channel loq is a pure state, hence 
5(a)f||a;o) < oo only if ujg = luq, i.e. £ cst — T. For- 
tunately, the identity channel is in accordance with the 
constraint |0) i— > |0), hence the estimation gives the iden- 
tity channel. In the second case it is straightforward 
to verify that the channel diag fullfils the constraints. 
Since iS(diag||£o) = is the minimal possible value we 
get £ est = diag. In the third case the support of ui n is 
a linear span of vectors = ( 1 00) + |ll))/\/2 and 

\4>+) = (|01) + |10))/\/2. Therefore, only for chan- 
nels with uj£ — a\<f>+){<f>+\ + b\ip + ){ip + \ the relative en- 
tropy is finite. However, the constraint requires that 
|0)(0| = (oT+M4)[|0)<0|] = a|0)(0| + 6|1)(1|, i.e. nec- 
essarily b = 0. In such case a = 1, because otherwise 
al is not a valid quantum channel. That is, only the 
identity channel satisfies the measured constraint, thus it 
minimizes the Kullback relative entropy. The estimation 
gives £ cst = I. In all these cases we find different esti- 
mations as in the unbiased maximum entropy approach 
(see Eq.(4.2)). The role of prior information in incom- 
plete process estimation deserves much deeper analysis 
than it is presented in these simple examples. However, 
such task is beyond the scope of this manuscript. 



Acknowledgments 

This work was supported by in part by the European 
Union projects QAP, by the Slovak Academy of Sciences 
via the project CE-PI and by the Slovak grant agency 
APVV and VEGA. 



7 



APPENDIX A: MAXENT SOLUTION FOR 

^proc = { h2 Q T ®a z ) £ } for qubit channels 

According to Section III this process observation level 
is equivalent to the following state observation level 

04 = {(2g T (E>a z } gi (a x (g>I) e ,(a y (g)I) e ,(a z (g>I) e } 
= {to, 0,0,0}. 

Maximum entropy estimation gives us the following state 
lo = ^cxp [-A • (ct® I) - 2dg T <g> a z ] (Al) 

where Z = Tr[cxp[— A ■ (a <g> I) — 2dg T ® er z ]] and A, d are 
Lagrange multipliers that can be determined by solving 
the system of algebraic equations 

= -InZ m = -—\nZ. (A2) 

8X dd 

Using the expression g T = ^(1 + Ft • <?) (Ft = 
(r x , —r y ,r z )) the state can be written in the form lo = 
\e~ R with 

R = A® |0)(0| + B® |1)(1| (A3) 
A = [(A + dr T ) -ff + dJ] 
B = [(X - df T ) ■ a - dl] . 



(<r x ) = to = the MaxEnt results in the state to = 

and, consequently, the estimated channel acts as follows 



1 



^e'-2Tr 1 [(C T ®/H = -7 



(A7) 



i.e., the whole Bloch sphere is transformed into the total 
mixture. 

For the case A = the first equation implies that either 
d = (leads to same solution as before) , or Ft = 0. If the 
the test state is chosen to be in total mixture (Ft = 0), 
the second equation leads to d = — arctanh(TO), hence 
the estimated state reads 



1 



LO 



4 cosh d 
1 



I ® (cosh d — sinh da z ) 



4 cosh d 
\l®\(I + ma z ). 



(A8) 



The corresponding process E est is given by the identity 
f[C] = dTr 1 [(e T ®/H,i.e., 



®T)(I®(I + mn z ))] (A9) 



T ri [e T ® (I + ma z )] 
(I + ma z ) . 



(A10) 
(All) 



Since the operators A <g> |0)(0| and B ® |1)(1| commute 
we can write 



e -R = e -A®|0)<0| e -.B®|l><l| 



(A4) 



= e 



|0)(0|+e- s ® |1)(1| 



Having in mind the operator identity 



= e x (cosh|^ + 



sinh 1 27*| 



-^—yv) 



\y\ 



(A5) 



we obtain 

Z = Tre- R = 2{e- d cosh |A + dr T \ + e d cosh \X - dr T \) . 
Inserting this expression into Eqs.(A2) we get 

a (A + dr T ) (A-dr T ) 



= S+ V U ' y + 5 1 - 

|A + drV| | A — drV 

(A + dr T ) ■ r T 



to(C+ + C_) = C+-C--2S. 



|A + dr T | 



<A6) 



where S± — e Td sinh |A ± dFr\ and C± — e^^coshlAi 
drT | • From the first of these equations it follows that A + 
dFr and A — dFr are collinear, i.e., A + drT — k(X — dFr)- 
This is possible only if either A = 0, or d = 0, or A = XFt- 
The case d = requires that A = 0, to = whatever 
test state Ft is used. Thus, measuring the mean value 



This transformation maps the whole state space into the 
single point £ = \ (I + ma z ). 

The last family of solutions of MaxEnt conditions is 
that the vectors A and Ft are collinear. In this case we 
reduced the number of unknown parameters to A = |A| 
and d. The first condition out of Eqs. (A6) then reads 



= e _ d sinh|(A + d)r| 



|A + d| 



(X + d) + e ^^^i (X-d) 



where we used r = \F\ = \Fr\- Analyzing all possible 
values for A ± d it follows that the absolute values can be 
omitted and the equations simplify to 



e d s. 



= e- d s+ 

e~ d c + - e d c 



TO = 



2re- d s+ 



e d c + + e d C- 



(A12) 
(A13) 



where s± = sinh[(A ± d)r], c± = cosh[(A ± d)r]. After a 
short algebra these equations can be rewritten into the 
form 



e Ar cosh[d(l 



er Xr cosh[d(l + r)] (A14) 



and 



m(e Xr cosh[d(l - r)] + e" Ar cosh[d(l +r)}) = 
= -e Xr sinh[d(l - r)] - e" Ar sinh[d(l + r)]) 
-2re- d sinh[(A + d)r]. (A15) 
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Unfortunately, we cannot give a general solution in a 
closed form. Consider therefore a special case and let us 
assume that the test state is pure, i.e., r = 1. In such 
case the solution reads 

1 1 1 — m 

d = -arctanh(— m) = - In ^—^ — (A16) 



A = - lncosh(2rf) . 



(A17) 



As a result we get 



uj = -e~ d [cosh(A + cf)/- sinh(A + d)f T ■ <?] <g> |0> <0| 
+ ^e d [cosh(A - d)I - sinh(A - d)r T ■ a] <g> |1)(1| 



In the language of Bloch vectors the transformation reads 



= (0,0,-m[l + f-r]). 



(A20) 



APPENDIX B: MAXENT SOLUTION FOR 

6 = {LI <8> 5%, La <8> I) e } = {™, 0} 

According to maximum entropy principle the 
state maximizing the entropy has the form 
u = ^ e -(\-3)®i-i®(fi-3) with 

Z = ^e^^' 3 ^ 1-1 ®^' 3 " 1 (Bl) 
= (Tre- t3 ){Tre- p - 3 ) = 4 cosh | A| cosh \fl\ . 



with Z = 2(e A + e~ A cosh(2d)). Let us denote by t the 
Bloch vector corresponding to a general input state £, 
then the estimated operation is given by the following 
prescription 

e^C' = 2Tn[(£ T ®/H 

= I(i + i m (i + r T -f T ))|0)(0| 

+ I(l-I m (l + f T .r T ))|l)(l| (A18) 

where we have used £ = |(1 +t ■ a). Taking into account 
that tr ■ tt — t ■ f we can write 

= \{I+\m{l+t-f)a z ). (A19) 



The values of A, jl are given by the following system of 
equations 

= -^^5 = -tanh|A|4- (B2) 
Z d\ |A| 



m = - — — m = -tanh /i — 
Z 9^ |^| 



and for the estimated state we get 



(B3) 



Z 4" 



^ = — ' ' ; ' = -I ®{I + m- a) . (B4) 



As a result we found that the estimated channel £ est maps 
the whole Bloch sphere into the point \ (I + m ■ a) . 
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